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ABSTRACT 

Researchers have reported (i) correlations of coronal mass ejection (CME) speeds and the total pho- 
tospheric magnetic flux swept out by flare ribbons in flare-associated eruptive events, and, separately, 
(ii) correlations of CME speeds and more rapid decay, with height, of magnetic fields in potential 
coronal models above eruption sites. Here, we compare the roles of both ribbon fluxes and the decay 
rates of overlying fields in a set of 16 eruptive events. We confirm previous results that higher CME 
speeds are associated with both larger ribbon fluxes and more rapidly decaying overlying fields. We 
find the association with ribbon fluxes to be weaker than a previous report, but stronger than the 
dependence on the decay rate of overlying fields. Since the photospheric ribbon flux is thought to 
approximate the amount of coronal magnetic flux reconnected during the event, the correlation of 
speeds with ribbon fluxes suggests that reconnection plays some role in accelerating CMEs. One pos¬ 
sibility is that reconnected fields that wrap around the rising ejection produce an increased outward 
hoop force, thereby increasing CME acceleration. The correlation of CME speeds with more rapidly 
decaying overlying fields might be caused by greater downward magnetic tension in stronger overlying 
fields, which could act as a source of drag on rising ejections. 
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1. INTRODUCTION 

In a coronal mass ejection (CME), magnetic forces in the low corona accelerate a hot (~ 1 MK) mass (^ 10^® g) of 
magnetized plasma at high speed (from a few hundred km sec“^ at the slow end to 2000 km sec“^ or more at the fast 
end) into interplanetary space. These events are the primary drivers of severe space weather disturbances at Earth 
(Gosling 1996), but key aspects of their initiation and subsequent evolution are not well understood. Prediction of GME 
speeds prior to their arrival at 1 AU is an important goal for space weather forecasts, both to accurately estimate their 
arrival times, and to account for the known correlation between GME speed and geoeffectiveness (Richciiusoii and Gane 
2011). 

Gharacterizing the forces that act upon CMEs as they begin to propagate away from the Sun is essential to un¬ 
derstand, and eventually predict, their dynamics. Here, we investigate two reported influences on the dynamics of 
flare-associated CMEs: (i) magnetic reconnection during the eruption process, as quantified by the amount of photo- 
spheric magnetic flux swept out by flare ribbons; and (ii) magnetic tension from large-scale magnetic fields that overlie 
the eruption site, parametrized by the rate of decay of horizontal field strengths with height, r, above the photosphere 
in potential field models. 

In the “standard model” of eruptive flares, dubbed the CSHKP (Carmichael-Sturrock-Hirayama-Kopp-Pneumann; 
e.g., Webb and Howard 2012; Svestka and Cliver 1992) model, magnetic reconnection in the corona beneath the rising 
ejection leads, either directly or indirectly, to particle acceleration. The resulting energetic flare particles then prop- 
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agate downward along the newly-reconnected magnetic fields toward the transition region and upper chromosphere, 
where they interact with the denser plasma and heat the atmosphere. This heating then powers dramatic intensity 
enhancements, which typically form two parallel, elongated emission structures, known as flare ribbons. These rib¬ 
bons always straddle polarity inversion lines (PILs), along which the radial photospheric magnetic field changes sign, 
which is consistent with ribbon emission occurring on conjugate footpoints of newly reconnected magnetic flux. In 
the simplistic CSHKP picture, the erupting structure itself originates directly above the PIL. Flare ribbons tend to 
systematically move apart as the ejection rises, as the reconnection links fields with footpoints that are farther and 
farther apart. Observers often study ribbons in Ha and UV images {e.g., in the 1600 A channel of the Transition 
Region and Coronal Explorer (TRACE) satellite (Handy et al. 1999). 

Flare ribbons overlie strongly magnetized areas of the photosphere within and near active regions (ARs), and it was 
long ago hypothesized that the rate of (unsigned) magnetic flux being swept out by the ribbons as they move apart 
should be related to the rate of coronal magnetic reconnection (Forbes and Priest 1984; Poletto and Kopp 1986). More 
recently, in a study of 13 CMEs with flares, Qiu and Yurchyshyii (2005) reported a very strong correlation between 
the total unsigned magnetic flux swept out by ribbons over the course of each event — which we hereafter refer to as 
the ribbon flux — and the CME speed: they found a linear correlation coefficient of 0.89. In a related study, Qiu et al. 
(2004) reported an association between the rate of flux being swept out by ribbons and CME acceleration. These 
findings are consistent: the time integral of acceleration yields velocity, and the time integral of the rate of magnetic 
flux being swept by ribbons (oc acceleration) yields the total magnetic flux swept by ribbons (oc velocity). 

The large-scale magnetic environment of CME source regions is also believed to influence CME dynamics. 
Kliem and Tdrok (2006) argued that CMEs are a manifestation of the torus instability, by which a toroidal cur¬ 
rent ring will be unstable to expansion if the external field surrounding the torus — which could inhibit any expansion 
via magnetic tension — decreases sufficiently rapidly in strength with radial distance from the center of the torus. 
Parametrizing the rate of decrease in external field strength with radius by defining a decay index, they theoretically 
characterized the dependence of the acceleration of the outward-expanding torus on the decay index, finding that 
a more rapid decay of external fields implies a more rapid expansion. Demoulin and Aulanier (2010) analytically 
modeled instabilities in more general configurations in terms of decay indices, and found that instability thresholds 
depended upon the held’s structure. Observationally, Wang and Zhang (2007) and Liu (2008) suggested that magnetic 
fields lying above flare sites can inhibit eruptions: using potential-field, source-surface (PFSS) models of magnetic 
fields above flare sites, they found that flares below horizontal fields that remained relatively strong with increasing 
height were less likely to be eruptive than flares below horizontal fields that decayed more rapidly with height. This 
is consistent with the idea that magnetic tension from stronger fields above flare sites is more capable of confining 
would-be ejections than weaker fields. Liu (2008) quantified the rate of decay of magnetic field with height using a 
decay index defined to be the fitted index of an assumed power-law decrease in horizontal magnetic field strength with 
radial height r above the flare site in the interval between about 45 Mm and 115 Mm above the Sun’s surface. 

Beyond the role of the rate of decay of overlying magnetic fields in the genesis of CMEs, overlying fields might 
also influence the speed of a CME after it has begun. When eruption occurs, it is physically plausible that magnetic 
tension from magnetic flux overlying a CME source region that is entrained by the leading edge of a rising ejection 
could slow the resulting CME. Also, flux above an eruption site that is not entrained would still need to be pushed 
aside, and might thereby impose an “effective drag” on the erupting CME, since this pushing would require doing work 
against magnetic pressure. Consistent with this picture, Xu et al. (2012) studied 38 CMEs and compared plane-of-sky 
CME speeds from the Large Angle Spectroscopic Coronagraph (LASCO; Brueckiier et al. 1995) CME catalog with 
decay indices derived from local. Green’s function potential field extrapolations above manually identified PILs. They 
reported that faster CMEs had larger indices; that is, field strengths above faster CMEs’ source regions decreased 
more rapidly with height. Although they did not compute correlation coefficients between their decay indices and 
CME speeds, their table contains the necessary information to do so, and we find linear and rank-order correlations 
of 0.64 and 0.60, respectively. We note that their extrapolations used evenly spaced grid points in height, but they 
apparently fitted the logarithm of height, so their fitted decay indices are probably biased toward the rate of change 
in field strengths at larger heights. To the extent that fields decay more rapidly at greater heights, as in their Figure 
1, this might bias their decay indices to be systematically larger than fits with points evenly spaced in the logarithm 
of height. 

We note that the approach of only modeling the field directly above the flaring PIL followed by Wang and Zhang 
(2007), Liu (2008), and Xu et al. (2012) does not account for the effects of complex magnetic structure in the erupting 
field. For instance. Sterling and Moore (2001) consider a scenario in which CMEs erupt sideways, from the edge of 
an active region. We also remark that the potential-field models that we analyze lack any electric currents, so must 


Reconnected Flux k. Overlying Fields vs. CME Speeds 


3 


certainly differ from the actual coronal fields, which evidently possess electric currents that store the energy needed to 
power CMEs. 

Based upon these reports of the roles of ribbon fluxes and decay indices in CME speeds, we set out to investigate 
the relative strengths of these two factors in determining CME speeds. Which factor matters most? And what is the 
unexplained variance in CME speeds after accounting for the joint influence of both factors? 

To address these and other questions, we analyze both ribbon fluxes and decay indices in a sample of 18 CMEs. 
In the next section, we present the events in our sample. Then, in Section -3, we describe our methods for analyzing 
each event, such as identifying ribbons and PILs, extrapolating fields, and determining decay indices. In Section 3.1, 
we illustrate these in a case study. We then present our results in Section 4, and conclude with a discussion their 
significance in Section o. 


2. DATA 

Once again, our aim is to quantitatively characterize the joint influence of ribbon flux and decay rate of overlying 
fields on CME speeds. To do so, we study a sample of 16 CME events from a seven year span, from 1998 to 2005. In 
Table 6, we show event times, solar source locations, CME speeds, and other parameters discussed below. 

In order to be able to compare our results with the findings of Qiu and Ynrchyshyn (2005), we analyze several events 
from that study. All the eruptive events we chose were “halo” CMEs, meaning that all CMEs were Earth-directed. 
Usually, halo CMEs originate from source active regions near disk center, so from the perspective of Earth-based 
observations, we can study them in greater detail (Webb and Howard 2012). More significantly, since halo CMEs are 
responsible for the strongest geomagnetic storms and other space weather phenomena (Webb and Howard 2012), it 
is very important for us to be able to predict them. Also, the sample we chose consists of mostly fast CMEs: only 
four have speeds less than 800 km s“^. We note that there are two events from the same source region, AR 10486, 
separated by several hours. These events are very famous, with remarkably high CME speeds (> 2000 km s“^) and 
ribbon fluxes. 

Ribbon fluxes and magnetic field maps in event source regions for all of these events were provided by Qiu and 
collaborators (Qiu and Yurchyshyn 2005; Qiu et al. 2007, and private communication). Full-disk, line-of-sight (LOS) 
magnetograms for these events were recorded by the Michelson Doppler Imager (MDI; Scherrer et al. 1995), and Qiu 
and collaborators extrapolated from the photospheric level to 2000 km, to estimate the flux density near the expected 
height of formation of ribbon emission. Here, the LOS field was used directly as an approximate measure of the 
radial magnetic field, and pixel areas were not corrected for foreshortening from the viewing angle. Images of flare 
ribbons were recorded either at Big Bear Solar Observatory, in Ha, or by the TRACE satellite, in UV. For each 
event, Qiu et al. generated 26 different ribbon masks {i.e., two-dimensional, binary maps) outlining the cumulative 
areas swept out by flare ribbons over the course of each event, by slightly varying selection criteria for ribbon pixels 
(Qiu and Yurchyshyn 2005; Qiu et al. 2007). These variations were undertaken to enable estimation of systematic 
errors due to mask selection parameters. Magnetogram data were interpolated onto the grid corresponding to the 
Ha/TRACE pixels. (This introduced artifacts at the edges of some magnetograms, which we discuss further below.) 
The lack of adjustment for foreshortening introduces errors of 10±7% over the set of all event fluxes; variations resulting 
from different ribbon-mask selection criteria are substantially larger. Note that, following Qiu and Yurchyshyn (2005), 
the ribbon flux is the total unsigned magnetic flux swept by the ribbon pixels — i.e., by including ribbon pixels in 
both polarities, the ribbon flux should be approximately double the amount of reconnected magnetic flux. 

3. MODELS AND METHODS 

Since we already have estimates of the ribbon flux and CME speed for each event, in order to compare the relative 
influence of decay indices on CME speeds we must still estimate the how the magnetic field changes with height above 
the CME source region for each event. 

As an initial step, we developed an automated (and therefore objective) method to define each CME’s source region 
on the photosphere from ribbon emission. It should be noted that CMEs are large-scale (possibly global) phenomena, 
so the small source region that we identify for each event should be interpreted more as the center or core of a much 
larger, surrounding source region than as an origin of all CME-associated dynamics. We expect that the magnetic 
field structure above the photospheric source region that we identify represents the large-scale magnetic environment 
of the erupting CME. 

As noted above, flare ribbons are thought to be the photospheric manifestation of coronal reconnection, so indicate 
which PIL(s) within the AR participated in the flare reconnection. Therefore, much as Xu et al. (2012) did, we used 
ribbon loci to determine, in part, where to model the structure of coronal fields. While ribbons typically straddle the 
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PIL, the eruption itself is assumed to originate above the PIL; hence, we define CME source areas to be near-PIL 
regions that are bracketed by flare ribbon emission. To identify near-PIL regions, we used an algorithm similar to 
that of Welsch and Li (2008), which requires both positive and negative fields with unsigned flux densities above a 
threshold of 100 Mx cm“^ to be close together: single-polarity magnetograms were created, then dilated by five pixels, 
and multiplied together; nonzero areas in the result correspond to the PIL mask. This flux-density threshold limits 
our analysis to strong-held PILs, by excluding insignihcant PILs from small-scale fluctuations in weak-held regions. 
We then found the intersection of these near-PIL regions with masks of ribbon emission. In Section 8.1, we study 
how varying parameters used in dehning near-PIL regions can effect our estimates of the rate of decay of overlying 
magnetic helds. 

In addition, for some of the events, (i) strong-held PILs and (ii) transient emission enhancements both occur far from 
the main haring PIL(s). While this spatially remote, transient emission is likely related to the hare, we do not expect 
the overlying magnetic helds above such distant regions to strongly inhuence dynamics of the erupting structure near 
the core ribbon areas. Here, we have implicitly assumed that the bulk of the CME originates from above the main 
haring PIL. We justify this assumption via Occam’s razor: hare ribbons are known to be associated with eruptions; we 
have no evidence that the eruptions that we study originated from parts of their source ARs remote from the largest 
area of ribbon emission; we therefore favor the simplest assumption given the observations, i.e., that the CMEs that 
we study originate in the corona above the main site of ribbon emission. Consequently, we choose to exlude helds 
overlying far-away, small-scale brightenings from our analysis. To do so, we kept only those pixels within a proximity 
threshold of 75 pixels from the geometric centroid of the ribbon mask for each event. Note that because pixel sizes vary 
between the interpolated magnetograms used in our extrapolations, this threshold does not represent a hxed distance, 
but was rather chosen to exclude transient emission that was clearly remote from the main hare ribbons. 

The resulting source area for each event is therefore the intersection of three sets: the ribbon mask; near-PIL regions; 
and pixels within the proximity threshold. We refer to pixels in this intersection as RPP pixels, for ribbon, PIL, and 
proximity. Since there are 26 masks for each event, there are in principle 26 sets of RPP pixels for each event. We 
note, however, that some choices of ribbon identihcation parameters lead to no overlap between our PIL masks and 
the ribbon masks. For this reason, there is only one RPP pixel in events ^6 and ^ 8 . 

We note that the decay index of the held above an individual pixel might not accurately represent the rate of decay 
of the held over the larger area through which an eruption must propagate. Consequently, increasing the number 
of pixels used to characterize the source region’s decay index might improve the correlation with CME speed. Set 
against this, however, including helds drawn from a larger area might introduce held structure that is irrelevant to 
the dynamics of some eruption(s), which could diminish the speed versus decay index correlation. For instance, while 
we have identihed all PIL regions, not all PILs participate in every eruption, so including the held structure above 
only those PILs near ribbons seems appropriate. As a practical matter, though, more pixels could be included easily 
enough, by separately dilating the ribbon and PIL masks (the “R” and hrst “P”, respectively, in “RPP”), which would 
likely increase the number of pixels in the intersection. In addition, the RPP intersection itself could be further dilated. 
In a case study, below, we do briehy investigate the sensitivity of the estimated decay index for one event to choices of 
masking parameters, but we leave a more thorough investigation of modihcations to our approach to a future study 

To model the helds overlying a CME’s source area, we extrapolate the potential magnetic held’s components above 
each RPP pixel. Like Xu et al. (2012), we use a Cartesian, Green’s function model. This approach treats the hux in each 
magnetogram pixel as a point source {e.g., Sakurai 1982), and assumes no sources exist outside of the magnetogram. 
The magnetograms are centered on the source ARs, and most are a few hundred seconds of arc on a side, meaning 
strong helds within several tens of megameters of RPP pixels are included in the extrapolations. 

No currents are present in our model held, but they are certainly present in the actual magnetic held in CME source 
regions. We expect, however, that current densities should usually decay rapidly with height, and, further, that the 
potential held approximates the actual held’s strength, if not its orientation. 

Next, we analyze the extrapolation(s) to characterize the overlying held’s structure. Following Liu (2008), we 
characterize the rate of decay of the overlying magnetic held with height by htting the extrapolated horizontal magnetic 
held strength, Rh = -|- R^, as a function of radial height r. For this ht, we assume Rh follows a power law in r, 

Rh oc , (1) 

where 7 is the decay index. We ht Rh(?’) for r G {46.5, 59.2, 72.3, 85.8,99.9,114.5} Mm. We have not investigated the 
orientation of the extrapolated held over this range, so the held morphology might not be arcade-like, as we assumed. 

This range of heights matches that htted by Liu (2008) in global, potential-held source-surface (PFSS) extrapolations. 
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These particular fitted heights, however, match those used by Suzuki, Welsch, and Li (2012) in fitting log(i ?0 + r), 
in which they are nearly uniformly spaced. Like Xu et al. (2012), however, we used Green’s function extrapolations 
instead of PFSS models. As noted above, Xu et al. (2012) fitted the logarithm of points linearly spaced in height r, 
which biased their fitted decay indices toward the decay rate at the large-r end of their fitting interval. Unfortunately, 
this bias from uneven sampling in log(r) did not occur to us until after much of our study had been completed. In 
fact, the heights we fit are not evenly spaced in either r or log(r); while they are more evenly spaced in log(r) than 
points evenly spaced in r would be, we still over sample the high-log(r) end of the range compared to the low-log(r) 
end. This can slightly weight our estimates of 7 by the rate of decay for larger r, but since our focus is on variations 
in 7 between events, a uniform bias should not be problematic for our analysis. 

We fit the horizontal magnetic field strength because we expect is a better measure of tension from the overlying 
field than the total field strength. A more negative value of 7 corresponds to a more rapid decline in as a function 
of height and thus a weaker overlying field. In our qualitative model, this would indicate that there is less downward 
magnetic tension due to the overlying field, and we would expect the CME to be faster than one from a source region 
with a less rapid fall-off in Rh- 

Suzuki. Welsch, and Li (2012) used two methods to find 7 , the fitted decay rate, from Bh{r). We use similar 
approaches here, estimating 7 two ways for each event and set of RPP pixels. In our Method 1, we first compute 
the average of transverse field strengths over the RPP pixels at each height, and then fit 7 from these average field 
values. In our Method 2, we first fit 7 from B\i(r) for each RPP pixel, and then average these 7 ’s to find a whole-RPP 
7 . Comparisons of 7 ’s computed two ways enables estimating the uncertainty in 7 . The results computed with each 
method are presented in Table 6 ; the significance of differences in these results will be discussed in Section . In 
both approaches, we used a linear, least-squares fit to the logarithms of both horizontal field strength, B^, and radial 
heights, r. For Method 1, we passed in the standard deviation of i?h over the RPP pixels at each height to the linear 
fitting procedure as an uncertainty weighting, to account for varying field strengths at that height. The uncertainty 
in 7 derived this way comes from the standard error of the estimation procedure. For Method 2, the uncertainty of 7 
comes from the standard deviation of the 7 values computed from all different pixels in the extrapolation region. 

In our conceptual model of the CME process (an extension of the CSHKP cartoon), reconnection of magnetic 
fields beneath the erupting structure might actively contribute to driving the CME’s eruption. This might occur 
either: directly, e.g., by momentum transfer to the CME from the “dipolarizing,” post-reconnection magnetic fields; or 
indirectly, e.g., by increasing the amount of flux that wraps around the ejection, thereby increasing the hoop force on 
the ejection. This reconnection can be measured by total photospheric flux swept out by flare ribbons. We calculated 
the magnetic flux by first taking the absolute values of the magnetograms and then summing over the product of these 
unsigned field strengths and the ribbon mask. Note that this double counts the reconnected magnetic flux; but since 
our analysis focuses on linear correlations, this factor of two is irrelevant. Since for each event we have 26 different 
ribbon masks, we used the average of the fluxes calculated from each of the masks and estimated the uncertainty of 
the flux using the standard deviation of these fluxes. The fluxes computed this way are listed in Table 6 . 

3.1. Case Study: A Typical Event 

In this section, we examine in detail the event on 26 July 2002 (^ 9), to illustrate how selection of RPP pixels 
and modeling of the overlying field was done, and to study the effect of different methods and parameters on our 
calculated value of 7 . The CME associated with this event had a speed of 818 km s“^ and ribbon flux 2.6 xlO^^ Mx. 
These are both around the average levels in our event list (although this speed is fast compared to most CMEs; e.g., 
Gopalswamy et al. 2010). We chose this event as a typical example of those in our list of CMEs to investigate the 
validity and robustness of our basic approach, as opposed to any particular features of this event. The magnetogram 
of this event exhibits some artifacts that are present in some other events’ magnetograms. We developed different 
extrapolation methods to test how different our results could be and how robust our original method is to these 
imperfections of our input data. 

The flare ribbon locations for this event, found from Ha flare observations, are shown (yellow contours) in Figure 
la, plotted over the corresponding magnetogram. In Figure lb, we show PIL regions (green contours) identified from 
the magnetogram. In Figure 2a, we show the 75-pixel proximity threshold (blue circle) centered on the mean position 
of the ribbon mask for this event. In Figure 2b, we show RPP pixels (red contours), which are the intersection of 
all three sets. Note that not all intersections of PIL regions and the ribbon mask are confined to a small area; some 
intersections are relatively far from the center of the ribbon emission, e.g., the relatively small patches of yellow at 
(a:, j/)=(275,215) and (325,120). As noted above, we do not expect the magnetic field directly above such small-scale, 
remote areas of ribbon emission to significantly influence the dynamics of the erupting CME, so they are excluded from 
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(b) PIL Mask 


Figure 1. (a) Contours of derived ribbon mask overlain on the interpolated full-disk MDI magnetogram for event ^ 9. (b) 
Contours of derived PIL regions over the same magnetogram. The magnetograms are displayed with a saturation level of 500 
Mx cm“^, and the white areas represent positive fields whereas black areas represent negative fields. 



Figure 2. (a) Contours of a circular mask of radius 75 pixels centered at the mean position of ribbon pixels, (b) RPP pixels 
(red) lie in the intersection of ribbon mask (yellow), PIL regions (green), and the circular mask (blue). The PIL region is 
selected with a threshold of 100 Mx cm~^. The potential field extrapolation was carried out above all red pixels. 


our analysis by the proximity threshold. Accordingly, we only include RPP pixels in the potential field extrapolation. 

Having defined the source region for this event, the extrapolation can be carried out. Using a Green’s function 
method, horizontal magnetic field strengths at points above the selected pixels were calculated from the line of sight 
field strengths of all other points in the magnetogram. This is done at heights of {46.5, 59.2, 72.3, 85.8, 99.9,114.5} Mm. 
According to our assumed power-law relation between CME speeds and the transverse field strengths, the power-law 
index 7 is —2.26 ± 0.2 using Method 1 and —2.09 ± 0.3 using Method 2. 

How robust is our estimate of this decay index? To find out, we conducted several tests, which we outline in the 
following paragraphs. 

To start, we tried some different parameters when deriving some of the key components for extrapolation to test the 
robustness of our methods. First, we note that when selecting the PILs, we used a rather high threshold 100 Mx cm“^ 
to exclude insignificant PILs. This substantially reduces the number of pixels identified for potential field calculation. 
To see if including more pixels would alter the result, we lowered the threshold to 50 Mx cm“^. Figure 3a shows the 
new PIL areas selected, which encompass many more pixels than with the higher threshold (Figure lb). Following the 
same extrapolation and fitting procedures, we determined the decay index from Method 1 to be -2.15, which is slightly 
lower than what we found with fewer pixels. This is expected since we now include some peripheral pixels with weaker 
fields. Still, this result is within one standard deviation of the original, indicating that original approach is robust. 

Also, we note that there are bad margins present in magnetograms for some of the events, including both this 
example and events #1, #3, #4, #5, #11, #12, #15 and #16. In these cases, field values are replicated from the 
edge of the estimated chromospheric field values to the margins of the magnetogram arrays. How significantly do these 
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Figure 3. (a) We show ribbon pixels (yellow), the proximity threshold (blue), PILs regions (green), and RPP pixels (red), 
but with the PIL regions determined using a different threshold: 50 Mx cm“^ instead of the 100 Mx cm“^ used in Figure 2b. 
Comparison with that figure shows that more RPP pixels are identified, including a new, small patch separated from the other 
patches, (b) The same quantities, but with the original PIL threshold, and the magnetic field artifacts at the magnetogram 
margins set to zero. Changes in neither the PIL threshold nor magnetogram margins significantly affect 7 . 


artifacts affect our estimated decay indices? To investigate the effect of bad margins on our decay index estimate for 
this particular event, we excluded the artifacts by setting the field there to zero (Figure -lb). The resulting 7 from 
Method 1 based on this modified magnetogram is -2.31 which is slightly higher than originally computed but also 
within one standard deviation. 

It is possible that the limited fields of view of our magnetograms could also affect our estimated decay indices, since 
fields outside of our magnetograms are ignored. But even within a magnetogram, how sensitively do our estimated 
decay indices depend upon magnetic flux relatively far from the flaring region? To see if our use of limited information 
matters, we decided to restrict our information to a circle with radius 150 pixels at the weighted event center. This 
specific radius was chosen because in the extrapolation for this event, the range of the heights are roughly comparable 
to 150 pixels and it is believed that, for the Green’s function extrapolation of a localized source distribution, the main 
contribution of the field at height H comes from fields with a circle of radius H centered at that point. A new test 
magnetogram was then obtained by setting the line-of-sight field outside the circle to zero. For this test, the decay 
index computed by Method 1 is -2.17, which is again only slightly different from our original value. This suggests 
that relatively distant magnetic flux should not significantly alter our estimates. These tests suggest that our original 
approach for estimating 7 values do not depend sensitively on our parameter choices or details of the magnetograms 
that we analyze. 


4. RESULTS 

Following the approach discussed above, we carried out potential field extrapolation for all 16 events, for all masks 
with nonzero ribbon pixels. We then computed the decay index based on these extrapolations. The speed of each 
CME, TCME, was taken from the CDAW CME list. A significant correlation between CME speeds (tcme), ribbon 
fluxes (i/j), and decay indices ( 7 ) was found. The CME speed data and ribbon flux yield a linear correlation coefficient 
of 0.76, with coefficient of determination 0.58; see Figure 5a. The scatter plot of tcme versus 7 in Figure ob also 
shows a clear correlation, with a linear correlation coefficient of -0.36, and coefficient of determination 0.13. A negative 
correlation between 7 — in which more negative values correspond to more rapidly decaying fields — and CME speeds 
accords with stronger overlying fields exerting stronger downward tension forces that lower the speeds of CMEs. To 
account for these mutual correlations, we compute the partial correlation between CME speed residuals and 7 , which 
excludes the effect of flux ip in both variables. This yields a correlation of -0.44 with coefficient of determination 0.19; 
see Figure 4. The existence of a residual correlation implies that extra information about CME speeds is contained in 
the decay index that is not accounted for by the ribbon flux. So, while it is plausible that 7 and ■i/' might be physically 
related (be., not physically independent variables), neither is redundant. 

Accordinly, we must characterize the joint dependence of ucme on 7 and '0. To do so, we fitted the CME speeds 
with a two-variable linear regression on ip and 7 , with all three variables standardized (i.e., the mean of each was 
removed, and the result was scaled by its standard deviation). The linear regression line was calculated to be 


i^CME = 0.73^ — 0.297 


(2) 
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Scatter Plot of Speed Residual vs. Gamma 



Gommo 

Figure 4. Scatter plot of CME speeds’ residual dependence on 7 , after removal of the dependence of both speeds and 7 upon 
ribbon finx 'ip. 


where wcme is the standardized speed and Tp is the standardized magnetic flux. (Note that our use of standardized 
variables here implies that no constant term is fitted in the model.) The multiple correlation coefficient — the 
correlation between the predicted and actual speeds — is 0.81, with a coefficient of determination of 0.66. For linearly 
related variables with Gaussian deviations, a coefficient of determination equal to one would imply that the dependent 
variable is completely determined by the independent variables; the difference between this coefficient and one is the 
fraction of variance unexplained by the independent variable(s). Under these assumptions, the increase in coefficient of 
determination from 0.58 for ijj alone to 0.66 for both ip and 7 represents a roughly 20% decrease in unexplained variance. 
Clearly, even with both variables, a substantial fraction of the variability in CME speeds remains unaccounted for. 

To characterize the sensitivity of our results to the particular set of 16 events that comprise our sample, we estimated 
the variances in our results using a delete-1 jackknife procedure (Efron 1982), i.e., by running our regression analysis 
on each of 16 subsets created from our event sample, each with one event deleted. We found the mean coefficient 
of pj in our regression to be 0.736 with standard deviation 0.034, and the mean coefficient of 7 to be -0.293 with 
standard deviation 0.062. Applying the delete-1 procedure to our standardized variables implies the means of each 
variable in the resulting 15-event subsets depart from zero, i.e., the subsets are no longer standardized. In this case, 
our regression procedure also generally returns a nonzero constant term. The mean and standard deviation of the 
fitted constant terms are 0.0036 and 0.045, respectively. Einally, we found the mean and standard deviation of the 
coefficient of determination to be 0.66 and 0.045, respectively. 

Our use of standardized variables implies that the independent variables in Equation (2) enter with equal weight, 
implying their predictive power can be directly compared from their regression coefficients. Consequently, the larger 
coefficient for ip, 0.73 versus 0.29 for 7 , implies that ribbon fluxes are much stronger predictors of CME speeds than 
decay indices. 

In Eigure 6 a, we plot the predicted CME speeds, based upon the regression in Equation (2), versus the actual CME 
speeds. The histogram of residuals between the actual CME speeds and the linear regression predictions is plotted in 
Eigure 6 b, computed in physical units. Discrepancies tend to be larger for very fast CMEs, compensated by a slight 
tendency for the model to overestimate speeds: there are 9 overestimates versus 7 underestimates. The outlier at the 
right end corresponds to event number 5. It is plausible that CME dynamics depend non-linearly on the variables we 
consider, instead of the linear dependence that we assume; but discriminating between the linear model and a more 
complex model would require a larger event sample. 

Also, inspired by the work of Robbrecht, Patsourakos, and Vourlidas (2009) on stealth CMEs —CMEs observed 
without “low coronal signatures” — we plotted flux and gamma as a function of CME speed to see if there were some 
lower bounds for flux and decay index below which there is no eruption at all. The flux versus CME speeds plot in 
Eigure oa suggests that there is no minimum ribbon flux required for CMEs, a result expected given reports of stealth 
CMEs. Usually such CMEs are slow (si 300 km s“^; Ma et al. 2010), and in order for a CME to be accelerated to 
a speed greater that the ambient speed of the solar wind (~ 450 km s“^), a significant amount of coronal magnetic 
reconnection, with consequent ribbon fluxes, appears necessary. In contrast, our fit of speeds versus 7 implies that, 
statistically, a minimum absolute value of 7 around 1.6 is required for CMEs. In our sample, though, several CMEs 
actually occur even with values of 7 less than this, highlighting that the fitted threshold is statistical. 
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(a) Scatter Plot of Flux vs. CME Speed 




Figure 5. (a) Ribbon fluxes as a function of CME speeds, (b) Gamma as a function of CME speeds. Event numbers are shown 
with each plotted symbol. 


5. DISCUSSION 

We have investigated the combined roles of (i) flare ribbon flux, ip, and (ii) the decay rate of overlying fields, 7 , 
in a set of 16 eruptive events. We have confirmed previous results that higher CME speeds are associated with both 
larger ribbon fluxes and more rapidly decaying overlying fields. Comparing the two effects, we found that CME speeds 
depend more strongly on ribbon fluxes than the rate of decay of overlying fields. 

Using a simple linear regression, we found that these two factors alone can explain 66 % of the variation in CME 
speeds. Given existing routine observations of (i) photospheric magnetic fields and (ii) EUV and/or Ha emission that 
can be used to identify flare ribbons, our results suggest that observers can readily make an approximate prediction 
of the speed of a flare-associated CME at the time of its eruption, prior to its appearance in a coronagraph. In 
terms of fractional errors in the predicted speeds of the relatively fast, halo CMEs in our sample, our mean \dv\/v 
is 25%, and our mean dv/v is -10%. Velocity estimates based on our results could be incorporated into existing 
CME propagation models {e.g., Zhao and Dryer 2014). Independent of variations in CME speeds due to acceleration 
occurring in interplanetary space, our errors in predicting CME speeds imply mean signed and unsigned errors in CME 
arrival times of 4 h and 10 h, respectively. This compares favorably with the current approach of fitting displacements 
of CME structure in LASCO image sequences {e.g., Mays et al. 2015), although we have only analyzed a small sample. 

We note that, compared to the study by Qiu and Yurchyshyn (2005), we included some additional events. In our 
study, we found that the correlation of CME speeds with ribbon fluxes was weaker than the correlation they reported. 
We hypothesize that either the effect of outliers in their smaller sample or new events within our larger sample are 
responsible for the difference. 

We remark that 7 values might depend significantly upon the method used to determine the potential held. Before 
we extrapolated helds using the Green’s function method, we hrst tried using PFSS models to estimate 7 above active 
regions in a sample of 47 CMEs. Following essentially the same PIL identihcation method discussed above (but without 
ribbon observations to identify which PIL regions were involved in these hares), we calculated the decay index 7 above 
PILs in these regions. In Table 1, we compare our results from the Green’s function method used described above 
with 7 ’s estimated from PFSS models. We hnd that Green’s function 7 ’s are systematically more negative than PFSS 
7 s, implying more rapid decrease of horizontal held strength with height. 

Given the similar approaches used, why do the inferred values for 7 differ systematically? First, it should be noted 
that these two sets of extrapolations were done above different sets pixels, which might by itself explain the differing 7 
values. But the difference might also be due to the lower spatial resolution of typical PFSS models (ours used spherical 
harmonics up to £ = 192), which yield smaller held strengths at low coronal heights. While our Green’s function 
method treats each pixel in the magnetogram as a single, delta-function source, PFSS models typically integrate over 
pixels to compute spherical harmonic coefficients. In this way, small bipoles (positive and negative hux close to each 
other) can get averaged out in low-resolution PFSS models, thus not contributing to the extrapolated held. With 
their better resolution, however. Green’s function methods can incorporate helds from small-scale bipoles into the 
extrapolations. Since the dipole moments from these small-scale helds are quite strong just above the photosphere 
but decay relatively rapidly with height, the held strengths computed using Green’s functions probably systematically 
decrease faster with height than PFSS helds, resulting in a more negative 7 . A more detailed study comparing methods 
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(o) Scatter Plot of Actual CME Speed vs. Predicted Speed 


Correlotion Coefficient = 0.81 


(b) Histogram of CME Speed Residual (Actual —Pred.) 
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Figure 6. (a) A scatter plot of actual CME speeds vs. predicted values, (b) Histogram of deviations of actual CME speeds 
from the linear regression predictions, in original nnits. 


of estimating decay indices would be worthwhile. 

Why do CME speeds depend more strongly on ribbon fluxes than decay indices? One possibility is that this is an 
artifact of our approach: since we only studied events with successful CME onset, we did not take into account the 
fact that there might be failed eruptions due to downward tension caused by overlying fields {e.g., Wang and Zhang 
2007; Liu 2008). In this sense, part of the effect of the decay rate of the overlying field on a CME eruption is obscured 
compared to the effect of ribbon fluxes: 7 might need to be negative enough for an eruption to be permitted at all, i.e., 
a high enough decay rate might be a necessary condition for an eruption in the first place. (Given that Xu et al. (2012) 
report eruptions for a very wide range of 7 values, however, it is not clear if any meaningful threshold in 7 exists.) 
Another possiblity is that potential field models are too inaccurate: perhaps discrepancies between the potential and 
actual magnetic fields above eruption sites are significant enough that our estimates for decay indices are wrong. Or 
perhaps decay indices themselves are irrelevant: maybe lateral motions in the early stages of eruptions enable CMEs 
to erupt “sideways” rather than pushing through overlying fields as we assume. 

Ultimately, it is clear that factors other than just the two that we consider also govern CME speeds. Likely 
possibilities include the internal magnetic structure of the pre-flare magnetic system that later erupts {e.g., the amount 
of flirx within this system, and / or the degree of its initial magnetic twist, perhaps quantified by magnetic helicity) 
and interplanetary magnetic field and density structures (including the effects of previous CMEs). We also remark 
that the CME speeds that we cite are plane-of-sky speeds; CMEs’ true speeds might be faster, which is an additional 
source of discrepancy between our simplistic model and the observations. 
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and The Catholic University of America in cooperation with the Naval Research Laboratory. SOHO is a project of 
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Table 1. Under CME Label, data sources are identified by H(Hq) and T(TRACE). We note that compared to 
Qiu and Yurchyshyn (2005) and Qiu et al. (2007), we add event and The information about date, time and speed 
of these events are from CD AW CME list, and the information about location and AR nnmber are from UH IfA active region 
label maps. TRACE’S 0.”5 pixels were summed 2x2. The * in the Decay indices (Method 2) of events #6 and #8 denotes 
that in the calcnlation of 7 for the two events, there is only one pixel in the RPP regions, thus precluding estimating standard 
deviations. 


CME 

Label 

Date 

Time 

[hh:mm:ss] 

Location 

Tome 
[ km s“^] 

AR 

# 

Gamma 

(Method 1) 

Gamma 

(Method 2) 

Flux 
[10^^ Mx] 

Pixel 

size (") 

1(H) 
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5(T) 
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6 (T) 
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10(T) 
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11(T) 
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12(T) 
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1.04 
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-2.05 ±0.5 

-2.09 ±0.2 
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1.04 


Table 2. In this table, we compare 7 ’s computed using Method 1 with our Green’s function method with 7 ’s from PESS models 
for an overlapping set of seven events. Apparently the magnitudes of 7 ’s computed from PESS models are either about the 
same or smaller that 7 ’s computed with Method 1. 


Event 7 , 

Number Green’s function 

I -1.32 ±0.09 
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